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1.0  INTRODUCTION 


Using  vertical  ionograms  to  quantitatively  determine 
radio  frequency  (RF)  propagation  conditions  is  a  mature  and  well 
established  technique.  In  principle,  a  large  number  of  scaled  vertical 
ionograms,  and  the  electron  density  profiles  that  can  be  '  drived  from 
them,  along  a  specific  path  of  interest  would  be  required  to 
completely  assess  the  RF  propagation  conditions  for  that  path. 
Typically,  the  midpoint  ionosphere  is  a  good  approximation  of  the 
ionosphere  along  a  specific  path  of  interest  and  RF  propagation 
conditions  for  that  path  can  be  determined  from  the  scaled  midpoint 
vertical  ionogram.  Previously  this  has  required  recording  a  vertical 
ionogram  at  the  desired  midpoint,  scaling  the  vertical  ionogram  and 
if  desired  inverting  the  scaled  vertical  ionogram  to  obtain  an  electron 
density  profile.  The  difficulty  of  placing  the  required  sounding 
hardware  at  the  midpoint  of  a  desired  path  can  be  resolved  by 
developing  a  system  capable  of  calculating  the  equivalent  scaled 
vertical  ionogram  at  the  midpoint  of  a  specific  path  of  interest  from 
any  oblique  ionogram  with  the  same  midpoint.  A  major  component 
of  the  Digital  Oblique  Remote  Ionospheric  Sensing  (DORIS)  program, 
is  the  development  of  an  automatic  scaling  algorithm  that  can 
produce  equivalent  scaled  midpoint  vertical  ionograms  from  oblique 
ionograms  produced  by  a  pair  of  Digisondes.  The  wide  range  of 
variations  that  can  reasonably  be  expected  to  be  present  in  oblique 
ionograms  recorded  during  disturbed  and  normal  ionospheric 
conditions  will  require  an  increased  level  of  sophistication  compared 
to  that  used  for  automatic  scaling  of  vertical  ionograms  recorded  at 
identical  ionospheric  conditions.  The  existence  of  a  reliable,  well 
tested  vertical  ionogram  scaling  algorithm  Automatic  Real-Time 
Ionogram  Scaling  and  True-Height  (ARTIST)  calculation,  as  part  of 
the  Digisonde  256  will,  however,  allow  adaptive  techniques  to  be 
utilized  in  the  development  of  the  oblique  ionogram  scaling 
algorithm.  While  the  specific  details  of  the  final  functional  oblique 
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ionogram  scaling  algorithm  will  depend  on  a  number  of  system 
factors  (e.g.  oblique  path  geography,  availability  of  polarization  data, 
transmit  power  and  waveforms),  application  of  adaptive  signal 
processing  and  pattern  recognition  techniques  will  provide  adequate 
flexibility  to  incorporate  the  present  off-line  oblique  ionogram 
scaling  algorithm  with  existing  Digisonde  256  software  and 
hardware. 
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2.0  APPROACH 


The  oblique  ionogram  scaling  algorithm  can  be 
represented  as  a  number  of  separate  but  interrelated  functional 
operations  including:  noise/interference  suppression;  trace 

enhancement;  trace  identification;  conversion  between  scaled  oblique 
echo  traces  and  the  equivalent  vertical  midpoint  echo  traces;  and 
echo  height  inversion.  The  interrelation  of  these  operations,  along 
with  the  required  input  data,  can  be  seen  in  Figure  2.1.  This 
algorithm  will  generate  scaled  oblique  ionograms  that  can  be  used  to 
produce  equivalent  scaled  midpoint  vertical  ionograms.  As  shown  in 
this  figure,  the  characteristics  of  the  ionosphere  at  one  or  both  ends 
of  the  path  of  interest  are  used  to  generate  an  estimate  of  the 
midpoint  ionosphere.  The  noise/interference  suppression  and  trace 
enhancement  operations  will  be  performed  via  a  two-stage  process 
utilizing  a  frequency  redundancy  technique  followed  by  an  adaptive 
Wiener  filter.  The  Wiener  filter  will  produce  a  collection  of  potential 
trace  segments,  from  the  received  oblique  ionogram.  The  trace 
identification  algorithm  will  then  determine  the  most  probable  set  of 
oblique  ionogram  traces  consistent  with  both  these  trace  segments 
and  additional  information  determined  from  vertical  ionograms 
obtained  at  the  receiving  node  of  the  oblique  path.  In  the  next  stage 
of  development  data  from  an  ionospheric  model  (at  present  we  are 
considering  the  International  Reference  Ionosphere  model),  will  be 
incorporated  with  the  vertical  ionogram  data  to  produce  an  improved 
estimate  of  the  midpoint  ionosphere. 


2.1  Frequency  Redundancy  Preprocessing 

The  first  noise/interference  suppression  operation 
performed  on  the  received  oblique  ionogram  utilizes  a  frequency 
redundancy  technique.  The  effectiveness  of  this  technique  is  based 
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on  the  fact  that  the  interfering  signals  produced  by  spurious  radio 
transmitters  are  manifested  as  large  amplitude  signals  at  every  slant 
range  for  a  narrow  frequency  band.  Generally  an  interfering 
transmitter  appears  at  only  one  sounding  frequency.  Since  the 
ionospheric  phenomena  of  interest  are  essentially  constant  over 
several  sounding  frequencies,  an  individual  sounding  frequency  (e.g. 
25  kHz)  step  with  the  highest  signal-to-noise  ratio  within  a 
predetermined  frequency  interval  (e.g.  100  kHz)  can  be  used  to 
represent  that  entire  interval.  The  result  is  a  relatively  noise/ 
interference-free  ionogram  with  information  content  equivalent  to 
an  ionogram  obtained  with  larger  (e.g.  100  kHz)  frequency  steps.  In 
the  preliminary  data  set  100  kHz  oblique  ionograms  were  generated 
from  50  kHz  oblique  sounding  data  by  selecting  the  50  kHz  channel 
with  the  highest  signal-to-noise  ratio  from  pairs  of  adjacent  50  kHz 
soundings.  In  the  proposed  operational  system  a  100  kHz  oblique 
ionogram  will  be  generated  from  25  kHz  sounding  data  by  selecting 
the  highest  signal-to-noise  ratio  channel  from  each  set  of  four 
adjacent  25  kHz  soundings. 


2,2  Noise/Interference  Suppression  and  Trace  Enhancement 

The  subsequent  noise/interference  suppression  and  trace 
enhancement  algorithm  uses  a  minimum-mean-square  Wiener 
filtering  technique  [Andrews  and  Hunt,  1969].  This  necessitated  the 
development  of  an  appropriate  mathematical  model  for  synthesis  of 
the  received  oblique  ionograms  from  the  leading  edge  components  of 
interest,  an  oblique-trace  spread  filter,  and  any  residual  noise/ 
interference  present  in  the  preprocessed  oblique  ionogram.  The 
details  of  this  model  are  illustrated  in  Figure  2.2.  The  preprocessed 
oblique  ionogram,  p0(t,f),  is  modeled  as  the  sum  of  a  noise/ 
interference-free  "ideal"  oblique  ionogram,  i0(t,f),  and  a  noise/ 
interference  signal,  n(t,f).  The  "ideal"  oblique  ionogram,  i.e.  the 
oblique  trace  echo,  is  modeled  as  the  response  of  a  linear  shift- 
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invariant  oblique-trace  spread  filter,  with  impulse  response  hs(t,f),  to 
the  desired  leading  edge  of  the  oblique  echo  trace,  le0(t,f).  To 
implement  a  minimum-mean-square  Wiener  filter  that  can  produce 
an  optimal  estimate  of  the  leading  edge  of  the  oblique  ionogram  from 
the  preprocessed  oblique  ionogram,  the  noise/interference  and  the 
leading  edge  signals  are  assumed  to  be  stationary  random  processes 
of  known  power  spectral  density.  In  addition  the  oblique-trace 
spread  filter  impulse  response  hs(t,f)  must  be  known.  The  function 
hs(t,f)  generates  the  pulse  width  of  the  transmitted  RF  pulses  and  is 
easily  constructed.  The  average  noise  function  n(t,f),  or  the 
corresponding  power  spectral  density  Pn(cot,cof),  can  be  obtained  by 
averaging  a  number  of  oblique  ionograms  that  do  not  contain  echo 
traces  (i.e.  "ionograms"  obtained  with  the  DORIS  receive  subsystem 
but  with  the  DORIS  transmitter  off).  With  these  data  the  two- 
dimensional  frequency  domain  optimal  noise/interference 
suppression  and  leading  edge  enhancement  filter  is: 

Hs^<ot,o)f)  Pijtot,(of) 

2 

Hs(cot,cof)|  P^CO^COj-j  +  PjcO^tOf) 

where  Pie((ot,cof)  is  the  power  spectral  density  of  the  leading  edge  of 
the  oblique  ionogram  which  of  course  is  not  known,  Pn(cot,cof)  is  the 
power  spectral  density  of  the  residual  noise/interference  signal  and 
Hs(cot,cof)  is  the  two-dimensional  Fourier  transform  of  the  oblique- 
trace  spread  filter  impulse  response  hs(t,f).  It  is  tempting  to  use  a 
generic  function  Pie(wt,cof)  to  represent  all  oblique  ionograms  but  this 
would  limit  the  effectiveness  of  this  technique,  given  that  the  filter 
would  only  produce  a  minimum-mean-square  estimate  of  the  leading 
edge  signal  relative  to  all  possible  oblique  ionograms.  To  overcome 
this  drawback  an  adaptive  approach  is  used.  The  specific  differences 
between  oblique  ionograms  can  be  attributed  to  three  major  factors: 
diurnal  variations;  geographical  orientation  of  the  specific  oblique 
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path  and  ionospheric  disturbances.  Since  all  three  of  these  factors 
would  be  manifested  in  the  vertical  ionograms  obtained  at  the  ends 
of  the  oblique  path,  the  scaled  vertical  ionogram  at  the  receiving  end 
of  the  oblique  link  is  used  to  estimate  the  leading  edge  of  the  oblique 
ionogram  for  that  specific  path  at  a  specific  time.  An  improved 
estimate  of  Pie(£Ot,cof)  could  be  obtained  by  using  the  vertical 
ionograms  from  both  terminal  points.  This  would  increase  the 
software  and  communication  requirements  since  the  IONHT  message 
from  the  transmitter  site  would  have  to  be  sent  to  the  receiver  site 
via  the  Air  Weather  Service  Automated  Weather  Network. 


3.0  EVALUATION  OF  OBLIQUE  TRACE  ENHANCEMENT  ALGORITHM 


The  performance  of  the  oblique  ionogram  trace 
enhancement  algorithm  was  evaluated  using  a  set  of  ionograms 
obtained  during  days  90-98  (March  30  -  April  7)  of  1988.  The 
oblique  path  was  approximately  1800  km  long  with  the  transmitting 
Digisonde  256  in  Goose  Bay,  Canada  and  the  corresponding  receiving 
system  at  Millstone  Hill  Observatory  in  Westford,  MA.  The  Goose  Bay 
station  used  a  horizontal  log  periodic  transmit  antenna,  while  the 
Millstone  Hill  station  used  the  standard  Digisonde  seven  loop  receive 
antenna  array.  The  protocol  used  to  collect  data  produced  sets  of 
three  ionograms,  consisting  of  a  vertical  ionogram,  an  oblique 
ionogram,  and  a  vertical  ionogram,  15  minutes  apart. 
Synchronization  of  the  two  stations,  which  is  required  for.  oblique 
ionograms,  was  performed  manually.  The  vertical  ionograms 
scanned  the  frequency  range  from  1  to  10  MHz  with  100  kHz 

sounding  frequency  increments,  while  the  oblique  ionograms 
scanned  over  a  frequency  range  of  4  to  18  MHz,  using  50  kHz 

sounding  frequency  increments.  The  results  presented  in  Figures  4 
through  14  were  generated  using  an  Apollo  660  workstation,  with 
the  required  algorithms  implemented  in  Fortran. 

A  representative  data  set  is  presented  in  Figure  3.1.  The 
upper  panel  shows  the  amplitudes,  as  determined  by  the  Digisonde 
256  MMM  technique,  of  a  vertical-oblique-vertical  ionogram  set, 
obtained  at  the  times  indicated.  The  pseudocolor  intensity  scale  used 
for  all  amplitude  displays  is  presented  at  the  extreme  left  of  Figure 
3.6,  ranging  from  black  for  the  smallest  amplitude  to  white  for  the 
largest  amplitude.  The  lower  panel  of  Figure  3.1  shows  the 

corresponding  polarizations,  with  the  vertical  ordinary  polarization 
displayed  in  red,  the  vertical  extraordinary  polarization  displayed  in 
green;  the  oblique  ordinary  polarization  is  displayed  in  yellow  and 

the  oblique  extraordinary  polarization  displayed  in  blue. 
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Figure  3.1  Representative  test  data.  Upper  series  consists  of  the  amplitude,  in  Digisonde 
256MMM  format,  of  a  vertical,  oblique  and  vertical  ionogram  set  obtained  at  the  times 
indicated.  The  lower  series  are  the  corresponding  polarizations. 
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The  oblique  ionograms  also  contained  local  vertical  echo 
traces  that  were  removed  before  noise/interference  suppression  or 
trace  enhancement  was  performed.  An  estimate  of  the  vertical 
echoes  in  the  oblique  ionogram  was  obtained  by  performing  a  logical 
NAND  operation  on  all  range  bins  tagged  as  part  of  either  the 

ordinary  or  extraordinary  traces  of  the  adjacent  vertical  ionograms. 

This  estimate  was  used  as  a  mask  to  remove  the  vertical  component 
in  the  oblique  ionogram.  Subsequently  any  removed  range  bind  was 
replaced  with  the  most  probable  amplitude  for  that  frequency.  The 
results  of  this  process  are  shown  in  Figure  3.2.  The  third  hop  vertical 
component  in  the  oblique  ionogram  was  not  present  in  the 
corresponding  vertical  ionograms  and  hence  was  not  suppressed. 

The  effects  of  the  frequency  redundancy  preprocessing 

can  be  seen  in  Figure  3.3.  Each  set  of  50  kHz  sounding  frequencies 
was  examined  and  the  most  interference  free  frequency  was  used  to 
produce  the  100  kHz  ionogram  seen  on  the  right  of  Figure  3.3. 

The  implementation  of  the  Wiener  filter  required  an 

estimate  of  the  power  spectral  density  of  the  interfering 
transmitters.  This  was  obtained  by  recording  50  interference-only 
ionograms  over  a  period  of  five  days  approximately  two  hours  apart. 
This  data  base  for  estimating  Pn(cot,co)  ls  appropriate  if  the  total 
number  of  interferers  is  approximately  constant  with  time.  The 
respective  power  spectral  densities  were  computed  and  averaged 
together  to  produce  the  data  presented  in  Figure  3.4.  The  strong 
components  along  the  cof  axis  indicates  that  the  interfering 
transmitters  produce  uniform  amplitude  signals  at  all  range  bins  at  a 
particular  sounding  frequency. 

The  power  spectral  density  of  the  leading  edge  of  the 
oblique  echo  traces  was  computed  by  first  synthesizing  an  oblique 
leading  edge  signal  from  the  ARTIST  scaled  vertical  data,  via  the 
secant  law  [Davies,  1969],  for  the  vertical  ionogram  immediately 
prior  to  the  oblique  ionogram  of  interest.  For  a  given  frequency  fv 


Figure  3.2  Representative  ionogram  set  in  upper  panel.  Lower  panel  is  the  corresponding 
oblique  ionogram  with  the  local  vertical  ionogram  suppressed. 


Figure  3.3  Original  50  KHz  oblique  ionogram  and  corresponding  1 00  KHz  ionogram  produced 
via  frequency  redundancy  technique. 
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reflected  vertically  at  a  virtual  height  h\  the  corresponding  oblique 
frequency  fo  reflected  at  the  same  virtual  height  is  fG  =  fv  d>0,  where 
00  is  the  obliquity  factor  given  as  O0  =  k  sec  0,  O  being  the  zenith 
take  off  angle.  The  factor  k  varies  from  1.0  to  1.1  for  distances  from 
0  to  3000  km  [Wilder,  1955].  Considering  curved  earth  geometry  the 
angle  <t>  is  given  by 

sin  (e/2) 

0= - ; - 

1  +  —  -  COS  (0/  2) 
a 

where  a  is  the  earth's  radius  (6370  km),  0  is  the  angle  at  the  center 
of  the  earth  subtended  by  the  endpoints  of  the  oblique  path,  and  h'  is 
the  vertical  reflection  height.  The  magnitude  squared  of  the  two- 
dimensional  Fourier  transform  of  the  synthesized  oblique  traces  was 
used  as  the  desired  power  spectral  density.  Figures  3.5  and  3.6  show 
these  results  for  the  representative  data  set.  In  the  scaled  vertical 
and  oblique  data  in  Figure  3.5,  the  E  trace  is  presented  in  yellow  and 
the  F  trace  in  white.  The  data  in  Figure  3.7  allow  comparison  of  the 
actual  vertical  ionogram  with  the  ARTIST  scaled  data,  and  the  actual 
oblique  ionogram  with  the  synthesized  scaled  data.  The  synthesized 
oblique  data,  superimposed  on  the  actual  oblique  data  in  the  right 
panel,  is  quite  similar  in  shape  but  shifted  up  in  frequency  slightly. 
As  expected,  since  the  scaled  vertical  data  is  derived  from  the 
vertical  ordinary  trace,  the  synthesized  oblique  data  more  closely 
matches  the  corresponding  oblique  ordinary  trace. 

The  Wiener  filter  behavior  is  seen  in  Figure  3.8.  The 
interfering  transmitters  produce  the  same  strong  signal  along  the  cof 
axis  that  was  evident  in  the  interference/noise  power  spectral 
density.  After  processing  by  the  filter  (shown  in  the  lower  two 
center  panels)  the  oblique  traces  (upper  right  panel)  are  enhanced 
(i.e.  the  interference  has  been  suppressed  and  the  leading  edge  of  the 
oblique  ionogram  is  more  apparent).  In  the  expanded  view  of  the 


Figure  3.4  Power  spectral  density  of  interfering  transmitters  and  additional  noise. 


Figure  3.5  ARTIST  scaled  vertical  ionogram  and  corresponding  scaled  oblique  ionogram 
synthesized  via  secant  law. 


Figure  3.6  Synthesized  scaled  oblique  ionogram,  magnitude,  phase  and  magnitude  squared 
(power  spectral  density)  of  corresponding  two-dimensional  Fourier  transform  (clockwise  from 
upper  left  panel). 


Figure  3.7  Vertical  ionogram  and  corresponding  ARTIST  scaled  ionogram  (upper  and  lower 
left  panel  respectively),  preprocessed  oblique  ionogram  and  corresponding  scaled  oblique 
ionogram  synthesized  via  secant  law  (upper  and  lower  center  panel  respectively),  synthesized 
scaled  oblique  ionogram  superimposed  on  actual  oblique  ionogram  (upper  right  panel). 


Figure  3.8  Preprocessed  oblique  ionogram  and  corresponding  two  dimension  Fourier 
transform  magnitude  and  phase  (from  top,  left  panel),  Wiener  filtered  version  of  preprocessed 
oblique  ionogram  and  corresponding  two  dimension  Fourier  transform  magnitude  and  phase 
(from  top,  right  panel),  Wiener  filtered  version  superimposed  on  preprocessed  oblique 
ionogram,  magnitude  and  phase  of  Wiener  filter  (from  top,  center  panel). 


Figure  3.9  Preprocessed  oblique  ionogram  and  corresponding  Wiener  filtered  version  (upper 
and  lower  panels  respectively). 
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original  oblique  data  and  the  resulting  enhanced  trace  segments, 
presented  in  Figure  3.9,  an  undesirable  artifact  is  present  in  the 
enhanced  trace  data,  primarily  at  those  sounding  frequencies  where 
the  original  oblique  data  was  corrupted  by  the  local  vertical  echoes. 
In  addition  to  the  residual  signals  resulting  from  the  incomplete 
suppression  of  the  local  vertical  echoes,  these  sounding  frequencies 
also  contain  multiple  hop  oblique  signals.  Since  neither  of  these 
components  are  accounted  for  in  the  Wiener  filter  the  resulting 
"impulse-like"  noise  is  produced.  To  a  lesser  extent  other  low  signal- 
to-noise  regions  of  the  oblique  ionogram  also  exhibit  this  processing 
artifact.  A  nonlinear  post-processing  linked-list  filter  [Fu,  1974]  was 
developed  to  reduce  this  artifact.  The  filter  rejected  as  noise  any 
group  of  contiguous  pixels;  i.e.  pixels  adjacent  in  either  frequency  or 
virtual  height,  of  less  than  a  predetermined  threshold.  In  the  data 
presented  in  Figures  3.10  and  3.11,  contiguous  groups  of  three  or 
more  pixels  were  considered  as  potential  oblique  trace  segments. 


Figure  3.11  Linked-list  post-processed  oblique  ionogram  superimposed  on  actual  oblique 
ionogram  and  linked-list  post-processed  oblique  ionogram  (upper  and  lower  panels 
respectively). 
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4.0  FUTURE  DEVELOPMENT 


At  present,  the  noise/interference  and  trace 
enhancement  components  of  the  endpoint  vertical  ionogram  based 
oblique  scaling  algorithm  have  been  implemented  and  tested  using  a 
limited  set  of  oblique  ionograms.  A  more  extensive  set  of  oblique 
ionograms  will  be  collected  both  to  further  test  these  algorithms  and 
to  develop  the  trace  identification  algorithm.  The  trace  identification 
algorithm  will  use  syntactic,  also  referred  to  as  structural,  pattern 
recognition  techniques  [Fu,  1974]  to  determine  the  most  probable  set 
of  oblique  ionogram  traces  consistent  with  both  these  trace  segments 
and  additional  information  determined  from  vertical  ionograms 
obtained  at  one  or  both  ends  of  the  oblique  path.  One  of  the 
potential  weaknesses  of  the  endpoint  vertical  ionogram  based 
oblique  scaling  algorithm,  as  proposed  in  Figure  2.1,  is  that  for  long 
paths,  in  particular  those  that  cross  the  terminator,  there  will  be 
significant  differences  between  the  endpoint  and  midpoint 
ionosphere.  This  can  be  addressed  by  using  adaptive  techniques,  as 
shown  in  Figure  4.1.  In  this  approach,  a  model-predicted  endpoint 
ionosphere  is  compared  with  the  measured  endpoint  ionosphere. 
Based  on  this  comparison,  the  model  parameters  are  adjusted  to 
produce  results  consistent  with  the  measured  data.  This  endpoint 
model  correction  is  then  used  to  correct  the  model  ionosphere  at  the 
midpoint.  The  synthesized  oblique  traces  produced  in  this  manner 
can  then  be  used  in  place  of  the  synthesized  oblique  traces  generated 
by  converting  scaled  endpoint  echo  traces  via  the  modified  secant 
law. 


Standard  Model 
Parameters 

O - 


Scaled  Vertical 
Endpoint  lonogram 
O - 


O  Oblique  Path 
Endpoint  Location 


Q  Time 


Oblique  Path 

Midpoint 

Location 


Electron  Density 


Model 

Parameters 
Consistent  with 
Endpoint 


Predicted 
Midpoint 
Electron  Density 


Synthesized  Oblique 
i.cho  Traces 

O - 


Ray  Tracing 
Algorithm 


Figure  4.1  Adaptive  Model-Based  Oblique  Echo  Trace  Algorithm 
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5.0  DORIS  MILESTONES 


The  milestone  charts  for  the  development  of  the  oblique 
ionogram  scaling  algorithm  and  the  overall  DORIS  project  are  found 
in  Figures  5.1  and  5.2  respectively.  At  present,  the  frequency 
redundancy  preprocessing  algorithm,  noise  suppression  and  trace 
enhancement  algorithm  and  post  processing  linked  list  algorithm 
have  been  implemented  and  tested  as  described  in  Section  4.0.  A 
more  extensive  oblique  ionogram  data  set  is  being  collected,  and 
after  the  initial  oblique  echo  trace  identification  algorithm  is 
developed,  this  data  set  will  be  used  to  study  the  performance  of  the 
endpoint  vertical  ionogram  based  oblique  scaling  algorithm.  This 
algorithm  will  produce  equivalent  digital  midpoint  vertical 
ionograms  of  virtual  height  vs.  sounding  frequency  with  2.5  km 
height  resolution  and  100  kHz  frequency  resolution.  These  results 
will  be  compared  with  actual  vertical  ionograms  at  the  midpoint  of 
the  oblique  path,  that  will  also  be  recorded  during  the  collection  of 
the  oblique  ionogram  data  set.  The  oblique  path  used  for  these 
experiments  will  have  the  transmitting  Digisonde  256  at  Goose  Bay, 
Labrador,  the  receiving  system  at  Wallops  Island,  Virginia  (or  vice 
versa),  and  data  from  the  East  Coast  OTH-B  Radar  quasi-vertical 
ionosonde  (Moscow  to  Columbia  Falls,  Maine)  to  determine  the 
corresponding  midpoint  vertical  ionograms. 
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Figure  5.1  Milestone  Chart  for  Oblique  lonogram  Scaling  Algorithm 
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